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METHOD AND APPARATUS FOR DETERMINING THE CEREBRAL 
STATE OF A PATIENT WITH FAST RESPONSE 

BACKGROUND OF THE INVENTION 
5 The present invention relates to a method and apparatus for 

determining the cerebral state of a patient. One application of the method and 
apparatus is determining the extent of a hypnotic state of the patient resulting, for 
example, from the administration of an anesthetic agent. That extent is often 
termed the "depth of anesthesia." In the method and apparatus of the present 

10 invention, changes in the cerebral state can be accurately and quickly determined. 

In a simplistic definition, anesthesia is an artificially induced state of 
partial or total loss of sensation or pain, i.e. analgesia. For most medical procedures 
the loss of sensation is accompanied by a loss of consciousness on the part of a 
patient so that the patient is amnestic and is not aware of the procedure. 

15 The "depth of anesthesia" generally describes the extent to which 

consciousness is lost following administration of an anesthetic agent. As the 
magnitude of anesthetization, or depth of anesthesia, increases, an anesthetized 
patient typically fails to successively respond to spoken commands, loses the eyelid 
reflex, loses other reflexes, xmdergoes depression of vital signs, and the like. 

20 While loss of consciousness (hypnosis, amnesia) and the loss of 

sensation (analgesia) are significant features of anesthesia, it should be noted that 
balanced high quality anesthesia must also consider muscle relaxation, suppression 
of the autonomous nervous system, and blockade of the neuro muscular jxmction. 
Sufficient muscle relaxation is required to ensure optimal operating conditions for 

25 the surgeon manipulating the patient's tissue. The autonomous nervous system, if 
not suppressed, causes the patient to respond to surgical activity with a shock 
reaction that effects heavily on hemodynamics and the endocrine system. To keep 
the patient completely motionless, the neuro muscular junctions transmitting orders 
from the brain to the muscles of the body need to be blocked so that the body of the 

30 patient becomes completely paralyzed. 
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While fhe need to determine the state of all five components of 
anesthesia is widely recognized, ascertaining and quantifying the state of hypnosis 
or depth of anesthesia in a reliable, accurate, and quick manner has been, and is, the 
subject of extensive attention. One reason for this is its importance. If the 
5 anesthesia is not sufficiently deep, the patient may maintain or gain consciousness 
during a surgery, or other medical procedure, resulting in an extremely traumatic 
experience for the patient, anesthesiologist, and surgeon. On the other hand, 
excessively deep anesthesia reflects an unnecessary consumption of anesthetic 
agents, most of which are expensive. Anesthesia that is too deep requires increased 

10 medical supervision during the surgery recovery process and prolongs the period 
required for the patient to become completely free of the effects of the anesthetic 
agent. A second reason for the continuing study and attention being given to 
monitoring the hypnotic condition of a patient arises because of its difficulty: that 
is, anesthetic agents alter the activity and state of the patient's brain and these 

15 chfiuiges are not always easy to detect. 

A measure of the depth of anesthesia that may be used for research 
purposes is foxmd in an Observer's Assessment of Alertness and Sedation or OAAS. 
The OAAS determines the level of consciousness or, conversely, the depth of 
sedation or anesthesia, based on a patient's response to extemal stimuli. One such 

20 assessment that classifies the deptti of anesthesia in six levels, is summarized by the 
table below. The transition firom consciousness to unconsciousness may be deemed 
to occur when the OAAS score changes from level 3 to level 2. Level zero 
corresponds to a state of deep anesthesia in which the patient shows no response to 
a very painful stimulus. 

25 
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OAAS 
Score 


Distinctive Characteristics 


5 


Patient replies readily to spoken commands, eyes open, awake. 


4 


Patient is sedated, but replies to spoken commands, mild ptosis. 


3 


Patient ceases to reply to loud commands, eye lid reflect present. 


2 


Patient does not reply to spoken commands, no eye lid reflex. 


1 


Patient does not react to TOP stimulation (50 mA) with movement. 


0 


Patient does not react to tetanic stimulation with movement. 



"Ptosis" is a drooping of the upper eyelids. "TOP stimulation" 
("train-of-four") is a very short, painful electrical (50 mA) stimulus applied to the 
ulnar nerve in the arm of the patient, repeated foxir times to evaluate the intensity of 
5 muscular contraction. In "tetanic stimulation" the electrical current (50 mA) is 
applied continuously for a period of time, such as 5 seconds. The ulnar nerve is the 
nerve which, when pinched, gives rise to the well known "crazy or funny bone" 
effect. 

While useful for research and other purposes, an OAAS scale 
10 provides only a limited number of scaling levels and is limited in practical use 
because of the attention required from the anesdiesiologist and the use of painful 
stimuli. 

It has long been known that the neurological activity of the brain is 
reflected in biopotentials available on the surface of the brain and on the scalp. 

15 Thus, efforts to quantify the extent of anesthesia induced h3^nosis have tumed to a 
study of these biopotentials. The biopotential electrical signals are usually obtained 
by a pair, or plurality of pairs, of electrodes placed on the patient's scalp at locations 
designated by a recognized protocol and a set, or a plurality of sets or channels, of 
electrical signals are obtained from the electrodes. These signals are amplified and 

20 filtered. The recorded signals comprise an electroencephalogram or EEG. 

Among the purposes of filtering is to remove electromyographic 
(EMG) signals from the EEG signal. EMG signals result from muscle activity of 
the patient and will appear in electroencephalographic electrodes applied to the 
forehead or scalp of the patient. They are usually considered artifacts with respect 
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to the EEG signals. Since EMG signals characteristically have most of their energy 
in a frequency range (40 - 300 Hz) which is different than that of the EEG, major 
portions of the EMG signals can be separated from the EEG signal. 

A typical EEG is shown in Fig.l. A macro characteristic of EEG 
5 signal patterns is the existence of broadly defined low frequency rhythms or waves 
occurring in certain frequency bands. Four such bands are recognized: Delta (0.5 - 
3.5 Hz), Theta (3.5 - 7.0 Hz), Alpha (7.0 - 13.0 Hz) and Beta (13.0 - 32.0 Hz). 
Alpha waves are foimd during periods of wakeftilness and may disappear entirely 
during sleep. The higher frequency Beta waves are recorded during periods of 

10 intense activation of the central nervous system. The lower frequency Theta and 
Delta waves reflect drowsiness and periods of deep sleep. 

By analogy to the depth of sleep, it can be said that the frequency of 
the EEG will decrease as the depth of anesthesia increases, while the magnitude of 
the signal initially often increases. However, this gross characterization is too 

15 imprecise and unreliable to use as an indication of such a critical medical aspect as 
the extent of hypnosis. Further, EEG signal changes during ginesthesia may not 
ftiUy correlate with changes in the hypnotic state of the patient. For example, it has 
been reported that in a 12-18 Hz frequency band, EEG activity initially increases as 
anesthetic agents are adntiinistered and only thereafter decreases as anesthesia 

20 deepens. 

The foregoing circumstance has led to the investigation and use of 
other techniques to study EEG waveforms to ascertain the underlying condition of 
the brain, including the depth of anesthesia to which a patient is subjected. It will 
be immediately appreciated from Fig. 1 that EEG signals are highly random in 
25 nature. Unlike other biopotential signals, such as those of an electrocardiogram 
(EGG), an EEG normally has no obvious repetitive pattems, flie morphology and 
timing of which can be conveniently compared and analyzed. Nor does the shape 
of the EEG waveform correlate well to specific underlying events in the brain. 
Hence, except for certain phenomena, such as epileptic seizures, which are readily 
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apparent from visual inspection of an EEG, the indication of other conditions in the 
brain in the EEG is much more subtle. 

Prefatory to the use of other techniques, the EEG signals are subjected 
to analog to digital signal conversion by sequentially sampling the magnitude of the 
5 analog EEG signals and converting same to a series of digital data values. The 
sampling is typically carried out at a rate of 100 Hz or greater. The digital signals 
are stored in the magnetic or other storage medium of a computer and then 
subjected to further processing to ascertain the underlying state of the brain. Such 
processing typically uses sets of sequential EEG signal samples or data points 

10 representing a finite block of time, commonly termed an "epoch." The analysis of 
the data is usually carried out on a moving average basis employing a given epoch 
and a certain number of backward epochs. 

Some of the techniques by which EEG signals can be analyzed in an 
effort to determine tilie depth of anesthesia are well described in Ira J. Rampil, A 

15 Primer for EEG Signal Processing in Anesthesia^ Vol. 89, Anesthesiology No. 4, 
pgs. 980 et seq-, October 1998. 

One such technique is to examine, in some meaningful way, how the 
voltage of an EEG signal changes over time. Such an analysis is termed a "time- 
domain analysis." Because of its generally random nature, an EEG signal is not a 

20 deterministic signal. This means that it is not possible to exactly predict future 
values of the EEG from past values in the manner that, for example, the shapes of 
past QRS complexes in an ECG signal can be used to predict future values for 
analytical and diagnostic purposes. However, certain statistical characteristics of 
random signals, such as an EEG, can be determined and used for analytic purposes. 

25 

Time-domain based EEG analysis methods have not proven greatly 
successful in clinical applications since the results do not behave in a completely 
consistent manner. However, such methods have been reported in the use of an 
electrical power parameter derived from the time-domain EEG signal voltage to 
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control administration of an anes&etic agent. Combinations of time-domain based 
statistic parameters have been used to analyze EEG data. Efforts have also been 
made to use the number of times the EEG signal crosses the zero voltage level in a 
given period to analyze EEG signal data. 
5 Time-domain based analysis is however useful in the study and 

quantification of burst suppression in the EEG. During deep sleep or anestiiesia, 
the EEG time-domain signal may develop a pattem of activity which is 
characterized by alternating periods or "bursts" of normal, or high frequency and 
amplitude, voltage signals and periods of low or no voltage, which periods are 
10 termed those of "suppression." The extent of this phenomenon can be expressed as 
a "burst suppression ratio (BSR)" which is a time domain EEG parameter 
describing the time the EEG voltage is in the suppressed state as a fraction of the 
sampling period. 

A second approach to analyzing EEG waveforms examines signal 
15 activity as a function of frequency, i.e. a "frequency-domain analysis." It has long 
been recognized that complex waveforms, such as EEG signals, can be 
decomposed, or transformed, into a plurality, or spectrum, of simple sine or cosine 
waves of various frequencies, amplitudes, and phases. Frequency-domain spectra 
can be obtained from sequential time-domain EEG signal data by a Fourier 
20 transform. Frequency-domain analysis analyzes the spectrum of frequency signals 
obtained from the transform to determine characteristics and features occurring in 
wave forms having the various frequencies of the spectrum. The results of an EEG 
frequency-domain analysis are typically graphically displayed as a power versus 
frequency histogram in which frequency is graphed on the abscissa and power is 
25 graphed on the ordinate. 

Further efforts to obtain useful information from 
electroencephalograms have employed higher order analyses, including the 
bispectrum and trispectrum. The bispectrum, which measures the correlation of 
phase between two different frequency components and quantifies the relationships 



-6- 



among the xmderlying sinusoidal components of the EEG, has received considerable 
attention. The bispectrum specifically quantifies the relationship between sinusoids 
at two primary firequencies fj and f2 and a modulation component at the frequency 
fl H- £2- A strong phase relationship between fi, f2 and fj + f2 creates a large 
5 bispectral value for frequency fi + f2. However, because the calculation must be 
performed using complex number arithmetic for several thousand fj, f2, and fi + f2 
frequency combinations, the computations to obtain bispectral information are 
rather arduous. 

For clinical use, it is desirable to simplify the results of EEG signal 
10 analysis of the foregoing, and other types, into a workable parameter that can be 
used by an anesthesiologist in a clinical setting when attending the patient. Ideally, 
what is desired is a simple, single parameter or index that quantifies the deptii of 
anesthesia on a consistent, continuous scale extending from frill alertness to 
maximally deep, but reversible, hypnosis. To be fully useful such a scale should 
15 maintain its consistency, notwithstanding the differing pharmacological effects of 
different anesthetic agents, as well as the differing physiologies of different 
patients. 

Various such parameters for relating EEG signal data to the hypnotic 
state of the patient are discussed in the literature. Several use frequency domain 

20 power spectral analysis. These parameters include peak power frequency (PPF), 
median power frequency (MPF), and spectral edge frequency (SEF). A peak power 
frequency (PPF) parameter uses the frequency in a spectrum at which occurs the 
highest power in the sampled data as an indication of the depth of anesthesia. The 
median power frequency (MPF) parameter, as its name implies, uses the frequency 

25 that bisects the spectrum. In the same fashion, the spectral edge frequency uses the 
highest frequency in tiie EEG signal. A modification of the latter is the SEF 95 
parameter, which is the frequency below which 95% of the power in the spectrum 
resides. 

To improve the consistency of an indicator of the hypnotic state or 
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depth of anesthesia, several parameters are often employed in combination. For 
example, the spectral edge frequency (SEF) parameter may be combined with the 
time-domain burst suppression ratio (BSR) parameter to improve the consistency 
and accuracy with which the depth of anesthesia can be indicated. 
S While parameters of the foregoing types can detect changes in the 

EEG caused by anesthetic agents and hence are useful in determining the depth of 
anesthesia, they suffer from an inability to be calibrated to behavioral endpoints and 
because of their sensitivity to the different EEG patterns induced by different 
anesthetic agents. 

10 More complex combinations of parameters are described in U.S. 

Patents 4,907,597; 5,010,891; 5,320,109; and 5,458,117 to Nassib Chamoun or 
Chamoun et al. and are employed in the anesthesia monitor product made and sold 
by the assignee of the patents. Aspect Medical Systems of Fnuningham, MA. The 
patents describe various combinations of a time-domain subparameter and 

15 frequency-domain subparameters, including a high order spectral subparameter, to 
form a single variable, termed the bispectral index (BIS), that correlates behavioral 
assessments of sedation and hypnosis over a range of anesthesia for several 
anesthetic agents. Because of this ability, the Aspect Medical Systems product has 
found clinical acceptance. 

20 The bispectral index, BIS, consists of the following three 

subcomponents: SyncFastSlow, BetaRatio, and Burst Suppression. The calculation 
of the subparameter SyncFastSlow utilizes bispectral analysis in the frequency- 
domain. The SyncFastSlow parameter corresponds to the logarithm of the ratio of 
the sum of all bispectral peaks in the frequency range 0.5-47 Hz divided by the 

25 sum in the range 40 - 47 Hz. The bispectral information in the SyncFastSlow 
subparameter does not, by itself, give sufficient information over the range of 
hypnosis thus requiring combination with the other subparameters. The BetaRatio 
subparameter gives the logarithm of the power ratio in the frequency ranges 30-47 
Hz £ind 1 1 - 20 Hz. It is a frequency-domain parameter that has been found to work 
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best in light sedation. As noted above, in very deep levels of anesthesia, EEG 
signal contains data samples in which the EEG activity is suppressed. The Burst 
Suppression Ratio obtained from a time-domain analysis of the EEG signal 
describes the relative content of burst and suppression in the signal. The Burst 
5 Suppression Ratio is operative in deep anesthesia in which the suppression occurs. 

The resulting bispectral index, BIS, is a combination of these three 
subparameters. The combining algorithm weights the difTerent subparameters 
according to their range of best performance. While flie details of flie algorithm are 
unpublished and proprietary, it is known that different subparameters or 

10 combination of subparameters are employed depending on the level of hypnosis or 
depth of anesthesia. For example, light sedation, it is necessary to use the bispectral 
SyncFastSlow subparameter in conjunction witii ttie BetaRatio subparameter in 
order produce reliable results. For deep anesthesia it is necessary to combine the 
bispectral subparameter SyncFastSlow with the Burst Suppression Ratio 

15 subparameter to produce reliable results. The algorithm appears circuitous in that 
in order to make the proper combination of subparameters required to accurately 
determine the depth of anesthesia, the algorithm must know what the level of 
anesthesia is which, in turn, requires the proper subparameter combination. 

Certain paradoxical behavior of the bispectral index (BIS) has been 

20 reported. See Detsch, et al. ''Increasing Isqflurane Concentration may cause 
Paradoxical Increases in the EEG bispectral index in Surgical Patients^\ Br. J. 
Anaesth. 84 (2000), pgs. 33-37. Because the index uses a plurality of 
subparameters and combinations thereof in different regions of hypnosis, this 
behavior may occur when the hypnotic level of a patient is at a boxmdary of the 

25 regions, for example, in the range between "surgical levels" and "deep hypnosis," 

Further, computation of the bispectral index (BIS) parameter requires 
averaging several epochs of EEG data. Thus, this index may be not sufficiently fast 
to detect changes in the state of a patient as is required in the clinical situation. See, 
Baker, et al. Electroencephalographic Indices Related to Hypnosis and Amnesia 
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During Propofol Anaesthesia for Cardioversion, Anaesthesia and Intensive Care, 
Vol. 28, No.4, 2000. Hence, the BIS index may indicate recovery several seconds 
after a patient has already opened his/her eyes. This can be a serious problem in the 
use of the BIS index. By knowing the depth of anesthesia, the anesthesiologist can 
5 more precisely control the amount of anesthetic agent administered to a patient. 
Often this results in a reduction in the amount of agent administered. However, the 
lessened amount of anesthetic agent increases the risk that the patient will awaken 
during surgery. It is therefore essential that an anesthesiologist knows immediately 
if a patient is approaching consciousness out of the hypnotic state. 

10 A different approach to the analysis of electroencephalographic 

signals is to attempt to quantify the complexity of the highly random EEG signal for 
use as an indication of the depth of anesthesia. This approach is based on the 
premise that neuronal systems, such as those of the brain, have been shown to 
exhibit a variety of non-linear behaviors so that measures based on die non-linear 

15 dynamics of the EEG signal should allow direct insight into the state of the 
imderlying brain activity. 

There are a mmiber of concepts and analytical techniques directed to 
the complex nature of random and unpredictable signals. One such concept is 
entropy. Entropy, as a physical concept, describes the state of disorder of a 

20 physical system. When used in signal analysis, entropy addresses and describes the 
complexity, unpredictability, or randomness characteristics of a signal. In a simple 
example, a signal in which sequential values are alternately of one fixed magnitude 
and then of another fixed magnitude has an entropy of zero, i.e. the signal is totally 
predictable. A signal in which sequential values are generated by a random number 

25 generator has greater complexity and a higher entropy. 

Applying the concept of entropy to the brain, the premise is that when 
a person is awake, flie mind is fiiU of activity and hence the state of the brain is 
more non-linear, complex, and noise like. Since EEG signals reflect the underlying 
state of brain activity, this is reflected in relatively more "randomness" or 
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"complexity" in the EEG signal data, or, conversely, in a low level of "order." As a 
person falls asleep or is anesthetized, the brain function begins to lessen and 
becomes more orderly and regular. As the activity state of the brain changes, this is 
reflected in the EEG signals by a relative lowering of the "randomness" or 
5 "complexity" of the EEG signal data, or conversely, increasing "order" in the signal 
data. When a person is awake, the EEG data signals will have higher entropy and 
when the person is asleep the EEG signal data will have a lower entropy. 

With respect to anesthesia, an increasing body of evidence shows that 
EEG signal data contains more "order", i.e. less "randomness", and lower entropy, 

10 at higher concentrations of an anesthetic agent, i.e. greater depth of anesthesia, than 
at lower concentrations. At a lower concentration of anesthetic agent, the EEG 
signal has higher entropy. This is due, presumably, to lesser levels of brain activity 
in the former state than in the latter state. See ^^Stochastic complexity measures for 
physiological signal analysis^^ by LA. Rezek and SJ. Roberts in IEEE Transactions 

15 on Biomedical Engineering, Vol. 4, No. 9, September 1998 describing entropy 
measurement to a cut off frequency of 25 Hz and Bruhn, et al. '^Approximate 
Entropy as an Electroencephalographic Measure of Anesthetic Drug Effect during 
Desflurane Anesthesia^', Anesthesiology, 92 (2000), pgs. 715-726 describing 
entropy measurement in a frequency range of 0.5 to 32 Hz. See also H. Viertio-Oja 

20 et al. "iVew method to determine depth of anesthesia from EEG measurement' in J. 
Clin. Monitoring and Comp. Vol. 16 (2000) pg. 16 which reports that tiie transition 
from consciousness to unconsciousness takes place at a universal critical value of 
entropy which is independent of the patient. 

The pertinence of the concept of entropy to the conscious and 

25 unconscious states of the brain is also supported in recent theoretical work (see 
Steyn-Ross et aL, Phys. Rev. E60 1999, pgs. 7229-7311) which applies 
thermodynamic theory to the study of the brain. This work points to the conclusion 
that when a patient imdergoing anesthetization passes from the conscious state to 
the vmconscious state, a thermodynamic phase transition of the neural system of the 
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brain takes place which is roughly analogous to the phase change occurring when 
water freezes into ice. During the process of freezing, an amount of entropy, 
proportional to the latent heat of the process is removed so that water and ice have 
different entropies. The conscious and unconscious states of the brain may 
5 therefore similarly be expected to have distinct, different values of entropy. The 
premise that loss of consciousness can be regarded as analogous to a 
thermodynamic phase transition, lends further support to the concept of entropy as a 
fundamental characteristic of the cerebral state of the brain and to the use of entropy 
in determining depth of anesthesia as employing a quantity reflecting the basic 

10 mechanisms of the brain rather than derived phenomena, such as power spectra, 
reflecting those mechanisms. 

In sum, the following can be said. First, certain forms of entropy 
have generally been foimd to behave consistently as a function of anesthetic depth. 
See Bruhn et al. and H, E. Viertio-Oja et al. ^^Entropy of EEG signal is a robust 

15 index for depth ofhypnosis^\ Anesthesiology 93 (2000) A, pg. 1369. This warrants 
consideration of entropy as a natural and robust choice to characterize levels of 
hypnosis. Also, because entropy correlates with depth of anesthesia at all levels of 
anesthesia, it avoids the need to combine various subparameters as in the bispectral 
index (BIS). Second, the transition from consciousness to imconsciousness takes 

20 place at a critical level of entropy which is independent of the patient. See Viertio- 
Oja et al. in J. Clin. Monitoring and Computing. Thirdly, and of particular practical 
significance, recovery of a patient toward consciousness from anesthesia can often 
be predicted by a rise of entropy toward the critical level. 

A number of techniques and associated algorithms are avedlable for 

25 quantifying signal complexity, including those based on entropy, as described in the 
Rezek and Roberts article in IEEE Transactions on Biomedical Engineering article. 
One such algorithm is that which produces spectral entropy for which the entropy 
values are computed in frequency space. Another algorithm provides approximate 
entropy which is derived from the Kolmogorov-Sinai entropy formula and 
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computed in Taken's embedding space. See Steven M. Pincus, Igor M. Gladstone, 
and Richard A. Ehrenkranz, "^4 regularity statistic for medical data analysis^\ J. 
Clin. Monitoring 7 (1991), pgs. 335-345. A program for computing approximate 
entropy is set out in the Bruhn et al., article in Anesthesiology. The spectral 
entropy and approximate entropy techniques have foimd use in analyzing the 
complexity of EEG signal data. 

Another technique for non-linear analysis of highly random signals is 
expressed in Lempel-Ziv complexity in which flie complexity of a string of data 
points is given by the number of bytes needed to make the shortest possible 
computer program which is able to generate the string. See Abraham Lempel and 
Jacob Ziv, "0« the complexity of finite sequences'\ IEEE Trans., IT-22 (1976) pgs. 
75-81. 

A still further approach fliat may be applied to EEG signal analysis is 
fractal spectrum analysis based on chaos theory. In fractal spectrum analysis, the 
EEG signal is divided into a harmonic component and a fractal component. The 
harmonic component includes the simple frequencies whereas the fractal 
component contains the part which is invariant under scaling in time. It has been 
found that the fractal exponent Beta which corresponds to the frequency power law 
l/f1^ increases consistently in the course of deepening anesthesia. 

BRIEF SUMMARY OF THE INVENTION 
An object of the present invention is to provide an improved method 
and apparatus for accurately determining the cerebral state of a patient, including 
the hypnotic or consciousness state of a patient and the depth of anesthesia that a 
patient is experiencing. 

A particular object of the present invention is to provide such a 
method/apparatus that can rapidly make such determinations, especially when a 
patient is emerging to the conscious state from xmconsciousness. 

The gist of the present invention is to combine an effective measure 
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of the cerebral state of a patient derived from EEG signal data, preferably a 
complexity measurement such as spectral entropy or approximate entropy, with a 
more rapidly obtainable measure derived from EMG signal data and to use the 
combination as a cerebral state indication. When used as an indication of the 
5 hypnotic state, or depth of anesthesia, of the patient, the measure derived from the 
EMG signal data enhances and confirms the determination of the hypnotic state 
made using the EEG signal data and renders ascertaining changes in the hypnotic 
state of the patient more rapid. This is of particular advantage in alerting an 
attending anesthesiologist to the possibility that an anesthetized patient may shortly 
10 regain consciousness so that the anesthesiologist can take timely, appropriate 
action. The measure derived from the EMG signal data may comprise spectral 
power data. 

Both the EEG and EMG signal data are typically obtained from the 
same set of electrodes applied, for example, to flie forehead of the patient. The 

15 EEG signal component dominates the lower frequencies (up to about 30 Hz) 
contained in the biopotentials existing in the electrodes and EMG signal component 
dominates the higher frequencies (about 50 Hz and above). 

The presence of EMG signal data shows that the patient is conscious 
and thus can provide a rapid indication of the conscioiis-unconscious state of the 

20 patient. Importantly, because of the higher frequency of flie EMG data signal, flie 
sampling time can be significantly shorter than that required for the lower 
frequency EEG signal data. This allows the EMG data to be computed more 
frequently so that the overall diagnostic indicator can quickly indicate changes in 
the state of the patient. 

25 In one embodiment of the invention, the EEG signal data and the 

EMG signal data are separately analyzed and thereafter combined into a diagnostic 
index or indicator. As noted above, because of the celerity with which changes in 
the anesthetic state of the patient can be determined from the EMG data, the overall 
index can quickly inform the anesthesiologist of changes in the state of the patient. 
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In another embodiment of the present invention, the spectral range of 
the complexity computations, i.e. entropy computations, is widened to extend into 
the EMG range. Thus, the spectral range over which the complexity computations 
are carried out to provide an indicator may extend from some lower frequency of, 
5 for example 0.5 to 7 Hz, up to a frequency above 32 Hz. To filter out power line 
interference, the spectral range may be divided into bands with the elimination of 
frequencies around 50, 60 Hz and 100, 120 Hz. For example, in an embodiment in 
which the spectral range extends to approximately 150 Hz, a lower frequency band 
(0.5 - 47 Hz) will contain mostly EEG signal data while two upper bands (63 - 97 

10 Hz and 123 - 147 Hz) will include primarily EMG activity. The use of a widened 
frequency range does not require a division of the spectrum into two segments or 
does the first embodiment because all components in the widened frequency range 
are treated in the same manner. And, any boundary within the spectral range would 
be artificial since the frequency bands for the EEG and EMG signal data are 

15 overlapping. The same analytical techniques are thus used for all levels of hypnosis 
from conscious state down to deep anesthesia, the paradoxical behavior found with 
indicators employing a plurality of subparameters and mles of combination for 
various levels of anesthesia is avoided. 

Further, the complexity measurement obtained in this second 

20 embodiment of the invention can be updated as often as is permitted by the higher 
frequencies of the EMG signal data in the widened spectral range of the complexity 
computation. This will provide a very current indication to the anesthesiologist of 
the depth of anesthesia of the patient. 

The indicator obtained from the signal complexity computation over 

25 the widened spectral range can be used in conjunction with a complexity 
measurement obtained only from the EEG portions of flie frequency spectrum to 
provide useftil information to the anesthesiologist regarding what portion of the 
indicator comes from cerebral activity and what portion comes from muscle 
activity. 
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Various other features, objects, and advantages of the invention will be 
made apparent from the following detailed description and the drawings. 

BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWING 
S Fig 1 shows an electroencephalogram. 

Figs. 2a and 2b are graphs showing values of entropy as compared to 
the conventional OAAS scale for a patient receiving an anesthetic agent; 

Figs. 3a and 3b are graphs showing values of entropy of a patient at 
sxirgical levels of anesthesia; 
10 Figs. 4a and 4b are graphs showing values of entropy as compared to 

the conventional OAAS scale for a patient emerging from anesthesia; 

Figs. 5a, 5b, and 5c are comparative graphic showings of various 
techniques for analyzing EEG signals; 

Fig. 6 is a flow chart showing one embodiment of the present 

15 invention; 

Figs. 7a, 7b, and 7c are graphs showing OAAS levels, EEG entropy, 
and EMG amplitude, respectively; 

Fig. 8 is a flow chart showing another embodiment of the present 

invention; 

20 Figs. 9a, 9b, 9c, and 9d are graphs showing OAAS levels and 

corresponding EEG and EMG based indicators of the depth of anesthesia; 

Fig. 10 is a graph showing combined EEG and EMG entropy values; 

and 

Fig. 1 1 shows apparatus for carr3dng out the present invention. 

25 

DETAILED DESCRIPTION OF THE INVENTION 
In a preferred embodiment of the invention, relevant information with 
respect to depth of anesthesia is extracted from the EEG signal data by computing a 
parameter that characterizes the amount of disorder or complexity in the signal. 
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Suitable matiiematical techniques include, for example, spectral entropy described 
in Rezek and Roberts, approximate Kolmogorov-Sinai entropy described in Pincus 
et al., and complexity described in Lempel-Ziv. The use of spectral entropy is 
deemed advantageous because of its superior computational simplicity and is 
5 described below. 

Computation of the spectral entropy of a signal according to Rezek 
and Roberts includes four steps. The first is the power spectrum calculation. The 
Fourier transform X(fi) of the signal x(ti) is computed by the fast Fourier transform 
technique (FFT). The power spectrum P(fi) is calculated by squaring the 
10 amplitudes of each element X(fi) of the Fourier transform: 

p(fO=x(fi)*x*(fi) (1) 

where X*(fi) is the complex conjugate of the Fourier component X(fi). 
15 The power spectrum is then normalized. The normalized power 

spectrum Pn(fi) is computed by setting a normalization constant Cjx so that the sum 
of the normalized power spectrum over the selected firequency region [fi,f2] is 
equal to one: 

20 t,Pnifi)-C„f^Pif,)=l (2) 

In the siramiation step, the unnormalized spectrsil entropy 
corresponding to the frequency range [fi,f2] is computed as a simi 

tPn(fi)^ogi-^) (3) 

25 

Thereafter, the entropy value is normalized. In order to normalize the 
entropy value to range between 1 (maximum disorder) and 0 (complete order), the 
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value is finally divided by ttie factor log (N[fi,f2]) where N[fi,f2] is equal to the 
total mraiber of firequency components in the range [fi,f2]: 

5 In the original work by Rezek and Roberts, the frequency range 

considered was the range below f2 = 25 Hz, as it is generally assumed that most of 
the EEG activity is confined to the firequency band below approximately 30 Hz, 

Fig. 2 shows, as a function of time, values of entropy as computed 
above as compared to an OAAS scale for a patient receiving an anesthetic agent. 

10 Fig. 2a indicates the OAAS level as determined by an anesthesiologist 

attending the patient. As noted above, at an OAAS level 5, the patient is fixUy 
awake whereas at the OAAS level 0 corresponds to a deep state of anesthesia in 
which the patient shows no response to tetanic stimulation. Horizontal line 10 
indicates a level at which transition from the conscious to unconscious state is 

15 deemed to take place, i.e. between OAAS level 3 and OAAS level 2. 

In the example shown in Fig. 2, the attending anesthesiologist 
considers the patient to have moved firom OAAS level 5 to OAAS level 4 at about 
three minutes. At about four minutes, the patient is deemed to have dropped to 
OAAS level 3. 

20 Thereafter, at about four and a half minutes, the patient is deemed to 

have lost consciousness as by failing to respond to verbal commands and the loss of 
the eyelid reflex. This is evidenced in the change firom level 3 to below level 2 and 
the crossing of horizontal line 10, 

Fig. 2b shows a value of entropy computed firom five seconds of data 

25 as graph 20 and a value of entropy computed as median values of twelve sequential 
five second epochs (sixty seconds) of data as graph 30, As can be seen firom Fig. 
2b, as the consciousness of the patient decreases firom the commencing of 
monitoring, both graphs 20 and 30 similarly decrease and cross horizontal line 40 
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which identifies the entropy level that characterizes flie transition from the 
conscious state to the unconscious state. 

In accordance with a protocol for the OAAS in use, the 
anesthesiologist commences the application of TOF stimulations to determine the 
depth of anesthesia on the OAAS scale. In the case shown in Fig. 2, the 
stimulations cause the patient to regain consciousness at about eight minutes. 

It will be seen from Fig. 2 that graphs 20 and 30 follow, and provide 
an accurate indication of, the state of consciousness of the patient, as presented on 
the OAAS scale. 

Figs. 3a and 3b show the values of entropy at surgical levels of 
anesthesia, i.e. when the OAAS scale is zero as shown in Fig. 3 a. Horizontal line 
40 in Fig. 3 is the same as horizontal line 40 in Fig. 1 and comprises the entropic 
value forming the borderline between the conscious and unconscious states. 

Figs. 4a and 4b shows a rapid recovery of a patient firom surgical 
levels of anesthesia to consciousness. The rise in the values of entropy, informs the 
anesthesiologist of the approaching recovery to the conscious state. 

While the invention has been described as using spectral entropy, Fig. 
5a, 5b, and 5c illustrate an example of anesthesia induction and emergence showing 
the suitability of the complexity measiirements of approximate entropy and Lempel- 
Ziv complexity as well as spectral entropy to determine the depth of anestiiesia. 
Measurements made using both shorter and longer samples of signal data are 
shown. 

Other techniques for analyzing the EEG signal data can also be used, 
if desired, such as higher order frequency domain analysis including the bispectrum 
and trispectrum, firequency domain power spectral analysis, and combinations of 
analytical quantities, such as the bispectral index (BIS). 

Measurement of electromyographic (EMG) activity contained in the 
biopotentials in the electrodes on the forehead, or other region of the scalp, of the 
patient can provide useful information concerning the state of an anesthetized 
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patient. As the level of anestiiesia approaches inadequacy, a painful stimulus causes 
a contraction of flie frontalis muscle (frowning) which can be detected as peaks in 
the EMG amplitude of the signal obtained from the electrodes applied to the 
forehead of the patient. EMG activity exists as long as the muscles are not 
5 paralyzed. This reaction can often be observed substantially before the pain 
eventually brings the patient to consciousness. EMG signal data can thus provide 
an early warning sign for the anesthesiologist to increase the level of anesthetics in 
order to prevent consciousness and awareness diiring surgery. 

Further, due to the high frequency range of the primary portion of the 

10 EMG activity above 40 Hz, a comparatively small time window is sufficient for 
computations using these frequency components, so that changes in the EMG 
activity can be detected substantially faster than changes in the EEG signal. 
Specifically, most of the component in the EEG signal resulting from brain activity is 
contained in the frequency range below 30 Hz. In order to obtain a good estimate of 

15 this activity by any mathematical procedure, the length of the signal used for 
computations has to be sufficient. In practice, the lowest frequency of the band sets 
the size for the signal length. For example, for an EEG signal band from 0.5 Hz to 32 
Hz, a signal sample 60 seconds long is required to obtain at least 30 cycles for each 
component. This sets a lower limit to the response time for assessing changes in 

20 patient brain activity from the EEG signal data. That is, the frequency with which all 
components of the EEG signal data indicator can be computed and ftilly updated is 
about once every 60 seconds. By contrast, an EMG signal data in a 63 - 97 Hz band 
requires only 0.5 seconds of data to obtain 30 cycles. The EMG signal data can thus 
be ftilly updated every half second. Because it can be so quickly updated, the EMG 

25 signal data can therefore provide an early warning sign for the anesthesiologist to 
increase, for example, the level of anesthetic agent administered to the patient in order 
to prevent awareness during surgery. 

Fig. 6 is a flow chart showing the steps for producing an improved 
diagnostic indication using EEG signal data and more rapidly indicative EMG 
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signal data in accordance with one embodiment of the present invention. In step 
100, the signal data corresponding to the biopotentials appearing in the electrodes 
placed on the scalp of the patient is obtained. In step 110, the signal data is 
subjected to spectral decomposition. This, may, for example, be carried out using 
5 Fourier analysis. 

The spectra are then divided into those representing the low frequency 
portions of the measured signal, for example, less than 30 - SO Hz and those 
representing the high frequency portion of the measured signal for example, those 
representing frequencies of 50 Hz and above. 

10 Thereafter, the EEG spectrum estimate is processed at step 120 to 

compute a parameter indicative of the state of activity of the brain. As noted above, 
it is presently deemed preferable to vise a computation of entropy for this purpose. 
However, other quantifications such as fractal spectrum analysis, Lempel-Ziv 
complexity, or bispectral or multispectral analyses, such as the bispectral index 

15 (BIS), can be used for this purpose. The result of this computation is the provision 
of an indication of the state of activity of the brain at step 122. 

As noted above, in order to obtain a good estimate of this activity by 
the mathematical analyses described herein, the length of the signal used for 
computations has to be sufficient. In practice, the lowest frequency of the band sets 

20 the lower limit for the signal length. In the case of the EEG indication generated at 
step 120, the lower limit for the signal is approximately 60 seconds. This means the 
indication can only be completely updated by repeating steps 100, 110, 120, and 
122 every 60 seconds and sets a lower limit to the response time of the EEG 
indication for assessing the patient's cerebral state. 

25 A power spectrum of the EMG signal is obtained in step 124, as by 

obtaining an amplitude spectrum and thereafter squaring the values of the amplitude 
spectrum to create a power spectrum. The EMG power spectrum provides an 
indication of EMG activity in step 126. 

Due to tiie high frequency range of the EMG activity, for example. 
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above 30 Hz, a comparatively small time window, for example 0.5 seconds, is 
sufficient to compute the EMG amplitude. This means that changes in the EMG 
activity can be detected and the indicator updated by repeating steps 100, 110, 124, 
and 126 substantially faster than changes in the EEG indicator, as shown 
5 graphically at steps 124, 124a, 124b, etc. 

In the example used above, the EMG indicator can be completely 
updated at a repetition rate of every 0.5 seconds. For simplicity in data processing, 
the EEG indication will typically also be recomputed every 0.5 seconds. However, 
since the EEG indicator requires 60 seconds of data, each computation 120, 120a, 

10 120b, etc. will use 59.5 seconds of old EEG signal data and 0.5 seconds of new 
EEG signal data. Thus, the changes in the cerebral state of the patient contained in 
the EEG signal data will be reflected only more slowly in the indication produced in 
steps 120 and 122 than the changes contained in the EMG indication. 

The EEG indicator and the EMG indicator are combined in a 

15 diagnostic indicator or index in step 128. 

The indicators produced in steps 120, 122, 124, 126 and 128 may be 
subjected to statistical treatment, such as averaging, if desired. 

The combined indication provided by the diagnostic index of step 128 
thus provides both reliable information of the activity state of the brain, such as the 

20 level of hypnosis or depth of anesthesia as direcfly foxmd in the EEG signal data, 
while full advantage can be taken of the rapidly obtainable information included in 
the EMG component of the signal which is a more indirect indication of the 
cerebral state of the patient but is particularly useful in alerting the anesthesiologist 
to the emergence of a patient from anesthesia. 

25 The components of the diagnostic indicator or index described above 

are shown in Fig. 7. An anesthetic agent is administered as a bolus at time zero. The 
patient enters imconsciousness, as shown by an OAAS level below line 10 in Fig. 7a, 
thereafter emerges for a short period of time, responsive to stimulation or the lack of 
further anesthetic agents, and is thereafter rendered unconscious. Fig. 7b shows the 
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entropy indication as obtained from steps 120 and 122 of Fig. 6. Fig. 7c shows the 
EMG amplitude obtained from steps 126 and 128 as a root-mean-squared simi over 
the EMG range of the Fourier spectrum. Data for five seconds are shown as the 
jagged lines. The smoother lines indicate one minute median filtered values. 
5 The graph of entropy as it relates to the hypnotic state of tiie patent 

resembles that of Figs. 2 and 4. With respect to the EMG activity, diiring the first two 
minutes following time zero, there is considerable EMG activity indicating that the 
patient is awake. Thereafter, the EMG activity decreases as the patient becomes 
unconscious. The unequivocal and intmiediate indication that the patient has regained 

10 consciousness at the ten minute point given by the EMG amplitude is clearly apparent 
fix>m Fig. 7c and will be reflected in the diagnostic indicator provided in step 128. 
This will advise the anestiiesiologist that the patient is emerging from the anesthesia. 

In another embodiment of the invention, incorporation of EMG signal 
data into a diagnostic indicator or index can be obtained by widening the frequency 

15 range of spectral entropy computations to one which extends from the EEG range 
into the EMG range thereby to include both EEG and EMG signal data. Although 
some small amount of EMG activity starts from frequencies of 1 Hz or lower, most 
EMG activity is traditionally quantified in a higher frequency range, such as the 
range from 40 Hz to 300 Hz (with notch filters at multiples of 50 Hz / 60 Hz in 

20 order to filter out the power line interference). However, there exists no clear 
frequency boimdary for the EEG and EMG data signals and an intermediate 
frequency range between 30 Hz and 50 Hz contains both EEG and EMG 
components overlapping each other. 

In providing an entropy based indicator obtained from a frequency 

25 range containing botii EEG and EMG spectra, it is to be appreciated that this 
represents a departure from the customary expression of, particularly, EMG signal 
data. That is, as noted above in connection with the first embodiment of the present 
invention, EMG signal characteristics are customarily expressed as a voltage 
amplitude, for example as a root mean squared spectral amplitude, and the 
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amplitude of the voltage will vary as a result of variations in the EMG signal data. 
By contrast, EEG entropy is a dimensionless quantity which describes the amount 
of disorder in the signal. When viewed from the entropic standpoint, the effect of 
EMG activity on the EEG signal is to create high frequency noise which increases 
5 the entropy of the combined signal . Entropy varies from 0 to 1, and the values are 
independent of the amplitude of the signal. The entropic expression of EEG signal 
data and the amplitude e>q)ression of the EMG signal data thus present a formal 
incompatibility which is overcome in the present invention by the use of entropy as 
a characteristics of both signal data. 

10 Further, this embodiment of the invention represents a departure from 

the previous entropic treatment of EEG signal data which has been limited to the 
frequency range in which EEG predominates as compared with EMG. These 
frequency ranges are, for example, frequencies of 25 Hz and below (see Rezek et 
al.) or 32 Hz or below (see Bruhn et al.). The present invention contemplates the 

15 use of frequencies in a range extending from some lower frequency, for example 
0.5 Hz, to a higher frequency which is in excess of 32 Hz. 

The second embodiment of the invention is explained using a 
frequency range of 0.5 Hz to about 150 Hz. It is deemed preferable to divide the 
extended frequency range into three-bands: a 0.5 — 47 Hz band; a 63 — 97 Hz 

20 band; and a 123 - 147 Hz band. The range is divided into the three bands at tiiese 
frequencies in order to avoid the power line harmonics at 50/100 Hz or 60/120 Hz 
depending on the frequency of the alternating current power mains. The lowest 
band contains most of the EEG components, while the two upper bands include 
primarily EMG activity. 

25 Fig. 8 is a flow chart showing the steps of producing an improved 

diagnostic indication or index using a widened frequency range for the computation 
of spectral entropy in accordance with the second embodiment of the present 
invention. In step 200, the signal data corresponding to flie biopotential signals 
appearing in the electrodes placed on the scalp of the patient is obtained. In step 
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210, the signal data is subjected to spectral decomposition, as by using Fourier 
transformation. As noted above, the spectral decomposition may be carried over a 
frequency range encompassing both the EEG signal data and the EMG signal data, 
for example, approximately 0.5 Hz to approximately 150 Hz. Also as noted above, 
5 the data is divided in spectral bands in order to omit frequencies at those of the 
power mains. 

At step 212, the lower frequency band of the EEG-EMG spectral 
range as well as the higher frequency bands are processed to compute a measure, 
such as spectral entropy, indicative of the complexity of the EEG and EMG signal 

10 data and state of activity of the brain and the frontal muscles. A combined indicator 
or index is provided at step 214. As described above in connection with Fig, 6, 
because the EMG signal data is available for complete updating more frequently 
than the EEG signal data, updating the combined indicator at the repetition rate by 
which the EMG signal data parameter can be updated provides a rapid indication to 

15 the anesthesiologist of any changes in the hypnotic state of the patient. 

Variations of the signal amplitude due to variations in the electrode 
contact properties etc. affect the EEG and EMG signal data in the same proportion 
to each other, so the combined EEG-EMG entropy indicator is not affected by the 
signal amplitude. 

20 Also, the approach taken in the second embodiment of the invention 

is consistent with the fact that the EEG and EMG frequency bands are overlapping. 
Thus no artificial boundary for EEG and EMG regions has to be defined since all 
frequency components are treated in the same way, i.e., a determination of the 
complexity properties. However, while the present invention features the use of a 

25 single algorithm for the range of conscious states of a patient, it will be appreciated 
the changes to the algorithm may be required, based on signal to noise 
consideration. 

A parameter containing only EEG entropic data from the low 
frequency band may be used as a separate indicator as shown in steps 216 and 218. 
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This can be used in connection witii the EEG-EMG entropic indication to allow the 
anesthesiologist to determine what portion of the EEG-EMG entropic indicator 
comes from brain activity and what portion comes from muscle activity. Thus 
while a clinical anesthesiologist will probably find the combined EEG-EMG 
5 entropic indicator to be highly usefril, a researcher may well find a comparison of 
EEG signal data entropy and the combined EEG-EMG entropy obtained over the 
wide frequency range to be of interest. 

Similarly, a parameter containing the EMG entropic data from the 
higher frequency bands may be used as a separate parameter by computing the 

10 entropy at step 220 and providing the indicator at step 222. The various diagnostic 
indicators are collectively shown at step 224 in Fig. 8. 

While both embodiments of the invention have been described above 
as utilizing Fourier transformation to decompose the spectrum of signals contained 
in the electrode biopotentials, other transformations, such as various wavelet 

15 transformations, may be employed. In order to obtain the corresponding indicators 
for pure EEG entropy, EEG-EMG entropy, and/or EMG entropy, the basis 
fimctions are divided into two classes, one including the basis fimctions for EEG 
activity, and the other including the basis ftmctions for EMG activity. After this 
classification, the indices for pure EEG entropy and EEG-EMG entropy are 

20 computed as described above. It is plausible that for better separation of the EEG 
and EMG components, a set of basis functions including wavelets emulating the 
shapes of the EMG spikes may extract the EMG component more effectively. 

Fig. 9 illustrates, as a function of time, the behavior of a pure EEG 
entropy parameter and the combined EEG-EMG entropy parameter along with the 

25 depth of anesthesia, as evaluated by an anesthesiologist on an OAAS scale. Also 
shown is EMG signal data as conventionally expressed as root mean squared 
spectral amplitude. The jagged lines show values computed from 5 seconds of data, 
while the smoother lines show one-minute median filtered values. An anesthetic 
agent is administered as a bolus at time zero. 



-26- 



During the first two minutes, during which time the patient is awake, 
there is a lot of EMG activity present as can be seen firom the graph of Fig. 9d 
showing EMG ampHtude. As a result of the high EMG activity, the combined 
EEG-EMG entropy shown in Fig. 9b also shows relatively high values, confirming 
5 that the patient is awake. At about the two minute point the patient loses 

consciousness as shown in Fig 9a when the OAAS score falls below line 10. 
Simultaneously, the EMG activity largely disappears as shown in Fig. 9d. The 
EEG-EMG entropy indicator and the EEG entropy indicator, produced at steps 218 
and 226, respectively, of Fig. 8 follow each other down below line 40 demarcating 
10 the transition to unconsciousness indicating deepening hypnosis. See Figs. 9b and 
9c. 

At about 6-7 nunutes both the entropy values of Figs. 9b and 9c show 
that anesthesia starts to lighten. This is because no further hypnotics have been 
administered. At about ten minutes, the patient wakes up and significant EMG 

15 activity again emerges. The pure EEG entropy indicator show in Fig. 9c predicts 
and indicates the recovery, but the EMG effect in the combined EEG-EMG entropy 
indictor shown in Fig 9b makes this niore apparent in its crossing of line 40. 

Thereafter following the ten minute point, patient is subsequently 
anesthetized again. The EMG activity gradually disappears as shown in Fig. 9d and 

20 EMG-EEG entropy indicator reduces to a pure EEG entropy indicator, as can be 
seen by a comparison of Figs. 9b and 9c. 

It is advantageous to scale the EEG-entropy parameter obtained in 
step 226 of Fig. 8 and the EEG-EMG entropy parameter obtained in step 218 in 
such a way that they coincide exactly when the EMG activity ceases completely, as 

25 this scaling allows for simultaneous graphical representation of the two pieces of 
information. This can be done in the following way. Take the firequency range 0.5 
to 32 Hz containing mainly EEG signal data to be firequency range Rj. Take the 
frequency range 32 to 147 Hz, with the power line frequencies removed, to be 
frequency range R2 for mainly EMG signal data. 
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Perform for the firequency ranges [Ri] = [EEG] and [Ri] + [R2] = 
[EEG] + [EMG] the steps 1-3 of the Rezek algorithm (Eqs. (l)-(3))to obtain 

sM = ZPn(Jinogi-^) (5) 



and 



Assume for a moment that Pn(fi) = 0 for all fj within the range [R2]- 
In this case, the normalization constant Cn[Rl + R2] = Cn[Rl] and, consequently, 

10 the normalized spectra Pn [Ri + R2] and Pn [Ri] (step 2) are equal. It follows that 
the unnormalized spectral entropies give by Eqs. (5) and (6) are equal. If tiiese 
entropy values are now normalized according to Eq. (4), they will no longer be 
equal, because the normalization factor log(N[Ri + R2]) is obviously larger than 
log(N[Ri]). This situation can be corrected by redefining the normalized entropy 

1 5 Sn[Ri] in the following way: 

S \r]_ ^Qg(M^i]) ^ S[R,] 

logiN[R, +R^l log{N[R, ]) logiN[R, 

The normalized entropy Sn[Ri+R2] can be defined as previously stated in Eq. (4): 

20 

According to these definitions, the EEG-EMG entropy varies from 0 
to 1, whereas the pure EEG entropy varies from 0 to log(N[Ri])/ log(N[Ri + R2]) 
< 1. The two entropy values coincide when there is no EMG activity so that P(fi) = 
0 for all fi within the range \R2]. When EMG activity is present, EEG-EMG 
25 entropy is larger than the pure EEG entropy. 



-28- 



In practice, a digital filter is commonly used in processing the signal 
data obtained firom the patient's biopotentials. Due to the characteristics of such a 
filter the computed entropy of a completely random signal, i.e. white noise, is 
usually slightly less than 1. For this reason, the entropies may be multiplied by a 
constant value so as to maintain the above described normalization. 

While the foregoing is described using two firequency ranges, Rj, R2, 
it will be appreciated that tiiese ideas can be straightforwardly generalized to the 
case in which the nxraiber of frequency ranges is larger than two. 

Fig. 10 shows the resulting normalized EEG entropy Sn[Ri] (thick 
curve) together with the normalized EEG-EMG entropy S]Sf[Ri + R2] (thin curve). 
When presented in this manner, the EEG entropy gives reliable information of the 
trend behavior of the patient's brain activity, while the EEG-EMG entropy responds 
rapidly to fast changes. 

Instead of the EEG entropy indicator and the EEG-EMG entropy 
indicator shown in Fig. 9, it is possible to compute the corresponding entropy 
indicator for EMG activity alone as shown in Fig 8, steps 220, 222, and to use this 
together with the pure EEG entropy indicator obtained in steps 216, 218. However, 
some care must be taken with respect to this approach. When EMG activity ceases 
due to relaxation of the muscles, some noise is left in the EMG range of the 
spectrum. The entropy of the noise may be relatively high and give a falsely high 
value for the level of EMG activity. Therefore, when EMG activity is considered 
separately and the concept of entropy is used for computations, a noise level should 
be established below which the EMG signal is considered to be zero. 

Apparatus for carrying out the present invention is shown in Fig.ll. 
Electrodes 300 are applied to the head of the patient in a desired manner. 
Preferably, at least some of the electrodes are appUed to the forehead of the patient. 
At least one pair and usually a plurality of pairs of electrodes are utilized. The 
biopotentials appearing in the electrodes are received in conductors 302 and are 
collected into patient cable 304. 
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Cable 304 connects conductors 302 to protection circuit 306 which is 
operative in the event the patient is subjected to electro-surgery or cardiac 
defibrillation. Electro-surgery employs alternating current at radio firequencies, 
typically between 300 and 3000 Hz to cut tissue and cauterize bleeding blood 
5 vessels. A defibrillator delivers a short current pulse to arrest arrhythmia in the 
heart muscle. Either of tiliese occurrences will significantly affect the signals in 
conductors 302 and, particularly, the EEG portion of the signals is usually rejected 
for fiirther use in determining tiie cerebral state of the patient. 

The output of protection circuit 306 is amplified by amplifier 308 and 

10 subjected to analog to digital conversion in analog/digital converter 310. Thereafter 
the signals are provided to bandpass filter at filter 312 that removes noise and line 
firequency harmonics firom the signals. The output firom bandpass filter 312 is 
connected to artifact detector 314. 

Artifact detector 314 detects artifacts arising firom electrocardiac 

15 activity, and other sources. The output of artifact detector 314 is connected to 
computational unit 316 which carries out the steps of the methods described above 
and shown in Figs. 6 and 8 and produces an output of the type shown in Figs. 3, 4, 
5, 7, 9, and 10 in display 318. Or, the information may be presented in display 318 
in numerical form. Display 318 may also display other physiological data, such as 

20 electrocardiographic data, breath rate, pulse, blood pressure, etc., obtained firom 
other monitors. 

Also, while artifact detector 314 is used to remove artifacts, the 
presence of artifacts can also be dealt with in the signal processing occurring in 
computational unit 316. For example, eye movements have been found to create 
25 simultaneous spikes in both the lower frequency EEG signal data as well as in the 
higher frequency EMG signal data. Sensing of the presence simultaneous spikes in 
both these frequency bands may be deemed to result from eye movements and 
particularly EEG signal data containing such artifacts can be eliminated from use in 
making the determination of the cerebral state of the patient. The same is also true 
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when excessive muscular activity and corresponding large EMG signal data is 
present. 

The invention has been described above in connection with cerebral 
states induced by the administration of an anesthetic agent. However, it will be 
S appreciated that the method and apparatus may be used in connection with other 
physiological conditions which are reflected in EEG and EMG signal data obtained 
from a patient and with drugs other than anesthetic agents. It is therefore recognized 
that other equivalents, alternatives, and modifications aside firom those expressly 
stated, are possible and within the scope of the appended claims. 
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